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A uS i'KACT 

A methodology for determining- the reliability 
of mechanical components is given. The necessary 
design data are pointed out. Complex fatigue re- 
search machines are described, which are gen- 
erating the required fatigue strength data in terms 
of cycles to failure and endurance strength. The 
test loading is that of alternating bending combined 
with steady torque. The design data obtained are 
presented. The design by reliability methodology 
and data are applied to the design of an actual al- 
ternator rotor inner shaft, and the results are 
compared with those obtained by conventional 
design procedure. 


SYMBOLS 


A 

average area, in. 2 * 

b 

positive integer 

c 

specimen radius, in. 

c’ 

power of ten exponent 

D 

specimen major diameter, in. 

d 

specimen groove diameter, in. 

e 

natural log base, 2.718 + 

I 

moment of inertia, in. 4 

J 

polar moment of inertia, in. 4 

f, kw, s, t 

stress concentration factor: b, 
bending; f, fatigue reduction; kw, 
keyway; s, torsional shear; t, 
part geometry 

k k 

coefficient of kurtosis 

k l 

life cycles factor 

ksk 

coefficient of skewness 

Ky 

coefficient of variation 


*Barred symbols are also used to represent 
statistically distributed variables 


k a,b, c, d, e, f modifying factors 
L^ut, y tensile load (see S), klb 

M bending moment on specimen, in. -lb 

Nf cycles to failure 


n 

R 

r 

r' 

r s 

®a, m 


®b, utj y 

s e 

Se 

s a, m 

T 

Tr 

W p 

Z 

A 

e 

v 

s 

( T 

Ox 


CO 


number of specimens 
reliability 

specimen groove radius, in. 
coefficient of correlation 
stress ratio, s a /s m 

graphical value material strength 
(see s), psi 

tensile strength, psi: b, breaking; 

ut, ultimate; y, yield 
corrected endurance limit, psi 
estimated endurance limit of the 
rotating beam specimen, psi 
material stress, kpsi: a, alterna- 
nating; m, mean 
torque on specimen, lb-in. 
test results: s, success; f, failure 
pan weight, lb 

standardized normal variable 
tolerance, in. 
angle, deg 

normal mean estimate 
difference statistic, psi 
normal standard deviation 
component of bending stress, psi: 
C, compresion; t, tension; a, 
alternating; m, mean 
component of radial stress, psi 

component of shear stress, psi: 
m, mean; a, alternating 
infinity 
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PRESENT MECHANICAL DESIGN techniques depend 
heavily on "modifying factors" to compensate for ■ 
limited knowledge of materials and loads. Some of 
these "modifying factors" compensate for imperfect 
knowledge concerning such tilings as surface condi- 
tion, size, temperature, stress concenti'ation and 
miscellaneous effects. In the best cases they are 
based on empirical data obtained over a long period 
of time, but in many other cases there is little real 
data and the selection is based on the designer's feel 
at the time. A final factor is often applied which 
generally reflects all the uncertainties involved and 
is probably influenced. as much by the seriousness 
of the consequences if the part should fail as by 
anything else. The resulting- design is probably 
considerably heavier than necessary. 

Light designs are important in aerospace tech- 
nology. Lighter designs might be possible if the 
designs could be based on the statistical matching 
of the available strength to the imposed stress. 

The NASA- Lewis Research Center in conjunction 
with The University of Arizona has undertaken the 
task of developing a mechanical design methodology 
of this type. The design methodology being devel- 
oped under this research is trying to reduce the 
number of these modifying factors by improved lab- 
oratory testing methods. It is dealing with im- 
portant design variables as statistical distributions 
rather than single valued estimates. It also makes 
provisions for designing an a priori reliability into 
a mechanical component. ' 

In order to work out this methodology the fol- 
lowing information is required; 

(1) Statistical strength distributions for a 
selected material. 

(2) Statistical stress variations in a loading 
configuration. 

(3) Statistical methods for relating stress to 
strength. 

This paper discusses the analytical, experimental i 
and application status accomplished during the first 
phase of this research. 

A typical rotating part design problem is given 
in this paper to illustrate the conventional design 
technique. The same problem is reworked using 
the reliability design methodology developed during 
Phase I of this research to see if weight savings 
can be accomplished. 


ANALYSIS 

DETERMINE TI1E FAILURE GOVERNING 
STRENGTH - In the engineering application of the 
design-by-reliability methodology, one of the 
most difficult tasks is to determine the actual 
strength and stress distributions involved for a 
specific material in a specific application. Two 
methods for determining distribution parameters 
for material strengths have been developed thus 
far; (1) Random variable functions; and (2) Di- 
rect experimentation. Much work remains to be 
done in this area. 

Random Variable Functions - The method of 
random variable functions as applied to material 
strengths is concerned with finding a statistical 
function which accurately represents the materi- 
al’s strength. Considerable research effort has 
been expended by Kececioglu (1)* and others to 
determine which distribution defines experimen- 
tal fatigue data best. Distributions which have 
been studied are the normal, lognormal and 
Weibull. 

A procedure to convert readily available ten- 
sile strength data, such as that given in Table A -3 
of Ref. (14, p. 600), to statistical functions which 
accurately represent the material's endurance 
strength would be very useful. This table gives 
yield strength, Sy, and ultimate tensile strength, 
Sut- In many cases where endurance strength 
distribution data for AISI steels are not available, 
use can be made of the common practice, based 
on the work of Lipson and Juvinall (2, p. 162), to 
obtain the mean estimate for endurance limit for 
steels as: 

Mg i = 0.50 S ut ; when S ut S 200 kpsi 

Mg, = 100 kpsi; when S ut s 200 kpsi 

The standard deviation lor the normal distribution 
can be obtained by making use of a conclusion 
based on Kececioglu's mu 1 others' work (3,4). 
These investigators have found that a standard 
deviation of 7 percent of the unmodified endur- 
ance limit can be u. c - r -'' ''' 'ittle or no test data is 

^Numbers in parentheses desi' v *e»le U»fer- 
cnces at end of papm 
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available. The final expression for S e becomes 
>S*. 0-07 S^). r 

"Marin's work (5, p. 127) can be extended to 
obtain a machined part's distributional endurance 
limit S e by considering each parameter in Eq. (1) 
below as a random variable. 

®e = ^a^b^c^kl^e^fSe (1) 

Products of random variables can be calculated 
by various methods (6,7). A later section gives 
an example calculation of this type for the case j 
of kjp kjj and kg unit normal modifiers, kg is 
included by a different method and Sg is based 
on experimental data rather than obtained from 
tables as explained previously. 

A second procedure for determining the dis- 
tribution of the fatigue strength of materials has 
been worked out using Monte Carlo Techniques. 
This procedure also requires knowledge of the 
constants k a through kf with an estimate of 
Sg from existing data. Digital computing tech- 
niques have been developed to calculate the first 
four moments of a statistic formed by the product 
of 1000 draws from the normal distribution of the 
known variables (8, 9, 10 and 11). Results ob- 
tained from the Monte Carlo computer calcula- 
tions compare within 5 percent of those obtained 
by the Algebra of Normal Functions Method when 
the coefficients of variation are small. There- 
fore, this latter method provides a simple ac- 
ceptable method for estimating strength distri- 
butions under these conditions. 

Direct Experimentation - hi some cases it 
may not be possible to estimate the distribution 
parameters for the fatigue strength of a material. 
This almost always is due to a lack of statisti- 
cally meaningful test data and/or engineering ex- 
perience with a particular material. A laboratory 
test procedure by which fatigue strength for dif- 
ferent materials could be obtained is desired. 

The first requirement of the test procedure 
would be to design suitable test specimens. Each 
test specimen must be carefully prepared and 
inspected to assure that it reproduces as many 
of the essential features of mo part application 
as possible. Many times it is necessary to scale 
various test specimen parameters making an 
effort to retain the same stresses, stress con- 


centration factors, stress ratios and material 
conditions. A statistically meaningful number of 
' these test specimens is also required. 

A machine to test these specimens under 
conditions similar to the part application is also 
needed. The fatigue testing machines used in 
this research effort are explained later. Groups 
of specimens are tested to failure at fixed alter- 
nating stress levels and stress ratios. Fig. 1 
shows the distributional S-N diagrams obtained 
from this testing. Five alternating stress levels 
and an endurance strength stress level are de- 
termined. The specimens are subjected to fixed 
magnitudes of reversed bending and steady torque 
to keep the stress ratio (alternating/mean), r s , at 
fixed values. Fig. 1(a) is for r g = oo and Fig. 1(b) 
is for r g = 0.825 for the cycle s-to -failure and 
for r s = 1. 03 for endurance. 

The cycle s-to -failure distributions for a par- 
ticular stress ratio can be converted to the ma- 
terial's strength distribution at any desired cy- 
cles of life. The cumulative histogram is formed 
to show the percentage of specimens failing in 
each stress cell for a fixed cycle-life value. 

From this cumulative histogram, the material 
strength distribution can be found by using sta- 
tistical methods. The resulting statistical dis- 
tribution is the material's strength for a fixed 
number of cycles of life and stress ratio. A num- 
ber of strength distributions may then be trans- 
lated onto the Goodman diagram, such as shown 
in Fig. 2, and the experimental failure governing 
strength distribution surface may thus be ob- 
tained. Fig. 2 shows a theoretical Goodman 
complex-fatigue strength surface. It seems that 
to establish an acceptable Goodman fatigue 
strength surface, at least four distributions 
(r s = 3,0, 0.70 and 0.30, at a specified life 

cycle, Kl = 10 n cycles) would be required. This 
coupled with the static ultimate tensile strength 
distribution, or at r g = 0, gives five sections on 
the contour to define the strength surface. The 
contours shown in Figs. 1 and 2 are for the two 
parameter normal distribution. A detailed task 
is being consider-ed to determine if other distri- 
bustions would fit this data better. 

DETERMINE THE FAILURE GOVERNING 
STRESS - The problem of determining the failure 
governing stress is concerned with determining 
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which failure governing theory applies best for 
the material, loads, dimensions, stress concen- 
trations and the like for the part in service. Ro- 
tating shafts can be analyzed using the von Mises 
Hencky failure criterion (8, 12). This criterion 
is applicable to ductile steels (8, pp. 152-154), 
and is theoretically valid for the elastic region 
only. It is recognized that the steel used in this 
research fits the ductile requirement in the un- 
notched form. The elastic region constrain is 
not met either, as the specimens are tested to 
fracture. The justifications for its use are that 
(1) correlation does exist between the results " 
predicted by it and the experimental results (13, 
pp. 38-41) and; (2) it is the only theoretically 
formulated criterion which explicitly incorporates 
all stress components involved. 

Fig. 3 shows the rotating shaft surface stress 
element for the two-dimensional case. The equa- 
tions for components of bending, radial and shear 
stresses are also given in this figure. Using the 
stress components as defined in the stress dia- 
gram, and according to the von Mises-Hencky 
criterion (14, p. 154), the failure governing al- 
ternating stress, s a , is given by 




■2 - - —2 , „—2 \l/2 

xa ~ ^xa^ya + °"ya + 3r xya J 


and the failure governing stress s m by 


m 


■( 


°xm °xm°ym + °ym + 3r xym 


1/2 


The test specimens studied in this research 
are explained later. For this specimen config- 
uration, the following conditions apply: 


'xa 


M_ 

G 


^M’ g M> 

(0.098 cf 3 , 0. 000244 d 3 ) 


( 2 ) 


s m ~ 


(3) 1//2 t 
vo ' xym 


1/2 

m i 

j 

c 


l-TS^fp, u-p) 

(0. 196 d 3 , 0.000488 cf 3 ) 


( 3 ) 


A detailed example of the use of similar equa- 
tions for a specific case is given later. 

The magnitudes of the corresponding failure 
governing strength (stress-to-failure) components 
can be obtained from the three-dimensional 
Goodman diagram by using graphical construction 
techniques, shown in Fig. 10. Fig. 10 shows a 
probabilistic, three-dimensional Goodman dia- 
gram. The top portion of the Goodman strength 
surface is obtained by drawing three straight 
lines to represent a distributional (normal) sur- 
face. The mean estimate loci are drawn as solid 
lines and the plus or minus three times the stand- 
ard deviation loci are drawn as dashed lines. 
These lines extend through experimental data 
points or from an ordinate of (p g g , erg ) to the 
junction of (Mgy, ersy) and a Sut)- A second 

stress plane is needed to define Sa and S m . 

This surface is defined by the probability density 
axis and the line passing through the origin sub- 
tending an angle of 02 with the abcissa. The 
angle 02 is given by the application stress ratio: 


J m 




= tan 


-1 




The intersection of this line with the mean 
Goodman line defines S a and 


BRIDGE THE GAP BY RELIABILITY THE- 
ORY - The reliability of a component can be de- 
termined from the basic concept that a no -failure 
probability exists when a given strength, S, value 
is not exceeded by stress, s. The probability 
that a stress of value s exists in the interval 
s^ - ds/2 to S]_ + ds/2 is equal to the area of 
the element ds, or to A^ in Fig. 4, or 


4 



1 


p( :; i - y - s - s i + t) = f(Bl)ds = Al 

The probability of strength exceeding s-^ is 
equal to the shaded area. A 2 , or 


P(S > s x ) F 



f(S)dS = A 2 


The probability of no failure, that is, the reli- 
ability, at s-^ is the product of these two prob- 
abilities, or 


dft = f(s- L )ds x 



f(S)dS 


The component reliability would then be all prob- 
abilities of strength being greater than all pos- 
sible values of stress (13,15) or 


■f-L 


/ f(s) 

f f(S)dS 

-OO 

J s 


R 

The reliability can also be written as 

.S 


ds (4) 


R 


- f m f 


f(s)ds 


dS 


( 5 ) 


Eqs. (4) and (5) can now be used to calculate 
the reliability of any component whose f(s) and 
f(S) are known. These equations carry limits of 
integration applicable to distributions defined 
over the interval from to +°°. For functions 
defined in different intervals these limits should 
be replaced by the lowest and highest values that 
can be used. Convolution integral techniques, 
transform methods, or a computer program uti- 
lizing Simpson's Rule (13, pp. 177-197) maybe 
employed to evaluate these equations and thereby 
determine reliability. 

If the density functions f(s) and f(S), rep- 
resenting the stress and strength distributions, 
respectively, are Gaussian or normal, as shown 
in Fig. 5, then they may be expressed as 


f(s) 


-l/2[(s-s)/cr £ 




f(S) = 



e -l/2[(S-S)/or s ] 


2 


Reliability is given by all probabilities that 
strength is in excess of stress or that S - s > 0. 
Using the designation £ = S - s, reliability is 
given by all of the probabilities that £ > 0. Let 
h(£) be defined as the difference distribution of 
f(S) and f(s), and as f(S) and f(s) are normally 
distributed, then h(£) is normally distributed 
also (16, pp. 215-216) and is expressed by 




-l/2[(£-£)/a ? f 


-y/2n 


where £ = S - s 


? 2 1/2 

and = (og + up 

Reliability would then be given by all prob- 
abilities of £ being a positive value, hence 




e -l/2[(£-£)/u £ ] 2 d£ 


The relationship between h(£) and the stand- 
ardized normal distribution can be utilized to 
evaluate the above integral. The transformation 
relating £ and the standardized variable z may 
be used which is 

z B LlI 
. °£ 

The new limits of the integrand are 
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for 


£ = 0 , 


z = LjlI = _X 


and for £ = +°°, z 


±gJL L = +QQ 


also d£ = arg dz 


j 

j 




If these conditions are substituted into the j j 
reliability equation, the following result is ob- 
tained: 



; Consequently, the reliability of a component is 
given by the area under the standardized normal ; 
! density function from the value of i 


' £ 1 

! z = - — to z = +« I 

! 

I i 

: The value of this area may be obtained from the ; 

! tables of areas under the standardized normal | 
density function, available in many references, j 

i 

! APPARATUS AND PROCEDURE l 


i The combined bending-torsion fatigue- j 
; reliability research machines are designed to j 
| simulate a shaft in service. The objective of the: 
immediate research program is to examine the j 
fatigue life of specimens made of SAE 4340 steel 
under combined loadings. The machine is de- 
| signed to subject the specimens to a predeter- 
! mined range of combined reverse bending and 
steady torque. A number of specimens are tested 
under this range of conditions to obtain the sta- ; 

! tistical distribution of the material's strength, j 
Each fatigue machine consists of a two- ' 

! section, rotating shaft with a test specimen j 
I locked in the center, as shown in Fig. 6. The 
horizontal shaft is coupled at each end to allow 
; for relatively free deflection when the specimen i 
1 is loaded. A horsepower, 1800 rpm motor 
; powers the shaft. The bending load is applied to 


the specimen by means of two voiws, one each on 
two bearings located symmetrically about the 
specimen on two commercial U)ol holders. 
the shaft, the yokes are connected by a horizontal 
link, which concentrates the load at a single ver- 
tical link in the center. The vertical link is then 
connected to either a long or a short loading lever 
arm. These loading arms make possible the ap- 
plication of a great range of bending stresses hi 
the specimen groove, by means of pan weights 
applied at the end of the loading arm. One lb of 
pan weight is equal to approximately 2000 psi 
bending stress in the groove. The torque is 
applied by means of a commercial Infinit-Indexer 
which is located on the back shaft of the machine. 
Table 1 summarizes the operational specifica- 
tions for the complex fatigue research machines. 
The machine is capable of producing, holding and 
transmitting to the rotating specimen steady 
torques of up to 5400 lb-in. and reversed bending 
moments of up to 3450 in. -lb. Specimens of 
diameters up to 1 in. can be tested in this ma- 
chine by changing the collets in the specimen 
holder. 

The SNAP-8 turbine shaft design specifica- 
tions were used as a guide to make the research 
data directly applicable to Aerospace problems. 
The American Society of Testing Materials Stand- 
ards were used as a guide for testing (17). Fig. 7 
shows the notched test specimen details. Bar 
stock pieces 6 in. long where carefully machined 
to the dimensions and finish specified in Fig. 7(a). 
Each specimen has a ground-in circumferential 
groove and keyway. The reverse bending moment 
M, torque T, bending stress a xa , and shear 
stress r X y m , are also defined in Fig. 7(a). 

Each specimen is installed in the test ma- 
chine and the instrumentation is checked for zero 
adjustment and calibration. Appropriate bending 
moment is applied to the specimen by putting 
weights on the load pan. Torque is applied to the 
specimen through the Infinit-Indexer with a suit- 
able wrench. The timing clock is set to zero and 
the machine is started. When the specimen frac- 
tures a micro switch stops the clock and the ma- 
chine; See Fig. 7(b) for the actual specimen be- 
fore and after fracture. 
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EXPERIMENTAL RESULTS 

: ]** 
■: A total of 396 specimens have been tested. 

Only 238 of these specimens yielded useful re- ; 
suits: 150 specimens yielded cycle s-to -failure 
data, 68 specimens yielded endurance strength 
data, and 20 specimens yielded static strength 
data. The remaining 158 specimens were not 
included in the reported test data because the 
machine test conditions varied over too wide a 
range during the test, or the specimen broke at j 
the keyway rather than in the groove. 

Tables 1 to 5 present the data obtained from 
the laboratory testing machines. Table 2 shows | 
the cycle s-to -failure data for stress ratios of ; 

«> and 0.825. Tables 3 and 5 show the notched 
and unnotched tensile strength data for a stress ! 
ratio of zero. Table 5 shows the endurance | 
strength data for stress ratios of « and 1.03. 
Table 6 shows a summary of the best fit equa- 
tions for the cycle s-to -failure data as given in 1 
Tables 2 and 4. 

DESIGN EXAMPLES ' 

A typical rotating part design problem is 
explained below. The inner shaft for an alter- 
nator rotor (see Fig. 8) is designed using mod- 
ifying factors to reduce the material's endurance 
limit. The shaft diameter following this design 
procedure came out to be 5/8 in. 

The inner shaft for an alternator rotor is 
redesigned using the probabilistic design method- 
ology. The shaft diameter following the distri- 
butional design procedure came out to be 1/2 in. 
This decrease in shaft diameter over the original 
requirement results in about a 38 percent de- 
crease in weight of material required to fabricate 
this part. None of the original design require- 
ments were reduced. These requirements were 
met by matching material to an application through 
the use of statistical methods. 

CONVENTIONAL SHAFT DESIGN - The inner 
shaft shown in Fig. 8 is being designed to support 
an alternator rotor. The shaft is made of AISI 
4340 steel, heat treated to 340/380 Brinell hard- 
ness with a machined finish. The rotor subjects 
the shaft to a radial bending load of 136.8 lb and 
a torque of 1000 lb-in. when delivering full load 
electrical power. The spline drive can exert a 


maximum misalignment radial force of 3 lb. The 
hot gas turbine drive will raise the shaft temper- 
ature to 500° F. The shaft rotates at 12,000 rpm 

o 

and is to have a life of greater than 10 cycles 
corresponding to a 0.999 component reliability. 
What should the shaft diameter be? 

There are various methods that can be used 
for machine design problems in which the shaft 
diameter is unknown. For aerospace applications 
life and reliability are very important consider- 
ations. Shigley (14, p. 482) has shown that under 
these conditions it is wise to use the von Mises- 
Hencky-Goodman method. This method requires 
knowledge of material strengths and application 
stresses. 

The example design problem specifies that 
this shaft is to be made of AISI 4340 steel, heat 
treated to 340/380 Brinell hardness. Table A-3 
(14, p. 600) gives the tensile properties for 1 in. 
round, drawn 1000° F condition, 377 Brinell 
hardness material as Sy = 160 kpsi and S^ = 

187 kpsi. 

Following Joseph Marin's suggestion (5, 
p. 127) the endurance limit will be reduced by 
modifying factors. These modifying factors 
attempt to account for the differences between the 
machined part's application environment and the 
endurance limit defined by laboratory testing. A 
machined part's endurance limit can be calculated 
by using Eq. (7) below: 

s e = kakbMVlkekfSe (7) 

The modifying factors will be explained later. In 
many cases, when endurance strength data for 
AISI steels are not readily available, use can be 
made of the common practice, based on the work 
of Lipson and Juvinall (2), to estimate endurance 
limits for steels as: 

S^ = 0. 50 S ut ; when S u ^ =£ 200 kpsi 

For AISI 4340 steel 

Sg = 0. 50x1, 87X10 5 psi =93.5 kpsi 

The fatigue strength modifying factors for this 
problem can be obtained from existing literature. 
The surface finish factor, k a , (14, Fig. 5-26, p. 167) 
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is 0,65 for 8^=187 kpsi. The size factor, k^, 

(18, p. 119) is 0.85 for specimens up to 2 in. 

The reliability factor, k c , (14, Eq. 5-14, p. 169) is 

k c - 1 - 0 . 08 D 1 - 0.248 = 0.752 

The standard deviation multiplication factor D 
is 3.1 for a survival rate of 99.9 percent. This 
factor begins to recognize that many design pa- 
rameters are not fixed, discrete integers but 
are more correctly represented by probability 
density functions of various kinds. The proba- 
bilistic design methodology explained later at- 
tempts to develop this concept still further. The 
temperature factor k^ is somewhat harder to 
estimate as only meager data exists. Curves of 
mean stress versus alternating stress for A 286 
alloy at 70° F and 1000° F are available (19, 

Fig. 2-36, p. 106). Each curve represents 
failure by rupture in 500 hr at the temperature 
indicated. Interpolating from the data given in 
this figure 1^ = 0.867. 

The inner shaft, as shown in Fig. 8(a), would 
fail at the shoulder fillet just to the left of the 
keyway should the rotor bind for some reason. 

The shaft is driven by the hot gas turbine through 
a spline coupling. At full electrical load the 
shaft carries 1000 lb-in. of torque. The direction 
of torque is given by the arrow in the load dia- 
gram. The bending moment diagram is given in 
Fig. 8(b). The shoulder just to the left of the 
keyway is the most dangerous section because of 
stress concentration. The clockwise bending 
moment here calculates out to be 126 in. -lb. 

The minimum theoretical stress concentration 
factors at this shoulder are (14, Fig. A-12-8,9, 
pp. 616-617) 

K =1.11 and K, =1.37 
b max D max 

for a shaft range of 1/2 to 5/8 in. diameter and 
the 1/8 in. shoulder radius. What to do when 
stress concentration is caused by three effects 


has not been invosf igaleil much. bill il seems rea- 
sonable that the combined effect should be no 
greater than the product of each factor. Therefore, 

K t = 1. llxl. 37 = U 52 

l comb max 

Notch sensitivity, q, (14. Fig. 5-27, p. 171) based 
on data (20, pp. 296,298) for a notch radius of 
0.125 in. and £5^= 187 kpsi is 0.90. The max- 
imum fatigue stress concentration factor K> „„ 
is calculated by 

K f =l+q(K t -l) = l+0. 90(1. 52-1)= 1.47 

max L max 

Finally, the minimum stress concentration mod- 
ifying factor k emin is given by 



The miscellaneous modifying factor, such as 
residual stress, directional characteristics, in- 
ternal defects, corrosion, plating are all taken 
to be unity for this example case. * 

Using Eq. (7), the modified endurance limit 
for this case would be 

S e = 0.65X0. 85X0. 752X0. 867x0. 680X1. 0X9. 35X10 4 
= 2. 28xl0 4 psi 

The bending moment M at the shoulder is 
126 in. -lb. The bending stress, cr at the 
shoulder is calculated, as follows: 


M 126X32 1.28X10 3 



C 


The torque T is constant over the shaft 
length at 1000 lb-in. The torsional stress 

' ■ I II I . I . ! » '- 

Some of the references cited in this report 
discuss methods for handling miscellaneous mod- 
ifying factors - Ref. (14) being especially useful. 
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r „ in the unknown diameter section of the 

1X1 

shaft is: 


T 1000X16 5.1X10 3 



neglecting the contribution of the alternating shear 
stress due to the vertical shear at the critical 
section. 

Analysis has been made of the various fail- 
ure theories suitable for defining stress compo- 
nents under these conditions (8, p. 37, 12, p. 479). 
The von Mises-Hencky failure theory has gained 
recognition as it seems to explain laboratory ob- 
servations well. The von Mises-Hencky stress 
components s a and s m are given by 


ordinate (F =0,2. 28X10 1 ) and abscissa (S ut = 
1.87x10 s , 0). The right boundary is defined by a 
point (Sy = 1. 60x10 s , 0) and an angle (tan -1 G = 

80. 5°) as the abscissa scale was choosen to be 
six times the ordinate scale to impi’ove accuracy. 

A third straight line is constructed to define S a 
and S m for this case. This line is defined by a 
point and an angle - (0,0) and (tan -1 0.865 = 40.8°). 
The intersection of Line 1 and Line 3 define the 
point S a = 12. 3xl0 3 psi and ^ = 85. OxlO 3 psi. 

Conventional design practice presently is 
using a margin of safety of at least 100 percent to 
insure a safe design. This requires that cr xa 
meet the following constraint: 

S a _ 12, 3x10 3 

0x3 ~ 2.0 2.0 


_ / 2 2 2 ' 
s m - \°xm “ <J xm cr ym + °ym + 3r xym i 


/ 2 2 2 ' 
3 a _ (°xa “ °xa°ya + °ya + 3r xya) 


For this problem, a reasonable model sets 
ff xm> ff ym> °ya and T xya all equal to zero . 
Therefore, 


s m = ^ T xym 


8.83X10 V 

j3 


s a " °"xa 


1.28X10^ 

n3 



A conventional Goodman diagram can be pre- 
pared from this data. Alternating and mean stress 
are plotted as ordinate and abscissa, respectively, 
in Fig. 9. Defining stress-to -failure as material 
strength makes it possible to use the Goodman 
boundary to represent the strength of this mate- 
rial as modified by environmental factors. The 
top portion of the Goodman boundary is obtained 
by drawing a straight line between two points - 


Substituting values gives 

1 :- -~ Xl0 - . ^ 6. 15X10 3 
d 3 

d 3 > 1 • 28x . 1 .P. 3 = 0. 209 
6.15X10 3 

or 

d s 0. 590 in. ! 

Based on this results use a shaft diameter of 
19/32 in. The results for the conventional design 
method are tabulated in Table 7 for purposes of : 
comparison with the probabilistic design method. 

PROBABILISTIC SHAFT DESIGN - The al- 
ternator shaft designed in the previous section 
shall be redesigned using the probabilistic meth- 
odology. The probabilistic design methodology 
being developed under Grant NGR 03-002-044 is 
based on Marin's technique using a modified form 
of the von Mises-Hencky-Goodman diagram. The 
primary change is to consider the necessary 
design parameters as probability density functions 
and the Goodman diagram as surfaces. Consid- 
erable laboratory testing is also being done to 
define material strength surfaces in this form. 
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The example design problem specified that 
the alternator shaft is to be made of AISI 4340 
steel, heat-treated to 340/380 Brinell hardness. 

The laboratory test program explained in the main 
text of this report gives test data for SAE 4340 
steel which is similar to AISI 4340, condition C4, 
per MIL-S-5000B, heat treated to Rockwell C 
35/40 per MIL-H-6875, with minimum tempering 
temperature of 1000° F, inspected according to 
MIL-I-6868. These tests showed that = 
{2.56x10 s ; 2.5X10 3 ) for notched specimens and 
Sut = (1.78xl0 5 ; 2.5X10 3 ) for unnotched specimens 
(see Tables 3 and 4). This mean estimate is 0. 96 
times that given in Table A-3 of (14, p. 600). The 
unnotched yield strength given in Table 5 is Sy = 

(1. 71X10 3 ; 3xl0 3 ) which agrees with the value 
given in Table A-3 of (14) fairly well (10.7 percent) 
and serves as a check on the test program validity. 

Since the laboratory test program is gen- 
erating data for endurance limits under combined 
conditions of reverse bending and steady torque 
for various stress ratios, not as many modifying 
factors are required to determine the parts ap- 
plication endurance limit. The test specimens 
were carefully machined to control surface finish 
so k fl has been included in the test program. 

These tests were run on round specimens with 
D = 0.735 in. , d = 0.4975 in. and r = 0.150 in. ; 
therefore, lq-, is already included. The relia- 
bility factor k c is introduced differently based 
on the overlap range of stress into strength as 
explained in the body of this report. All of the 
test specimens run on this project so far have 
been in a laboratory ambient temperature of 75° 

±5° F; therefore, an estimate of the modifying 
factor for temperature k^ is required. One 
value for k^ = 0.867 can be obtained from (19), 

Fig. 2-36, p. 106 for unnotched A 286 alloy steel. 

A second value for k^ = 0. 921 can be obtained 
from (8), Fig. 1.20, p. 46 based on (20) for un- 
notched, 1.8 percent nickel, 0.8 percent chrom- 
ium alloy steel. To obtain a normal probability 
function, define the range of k^ as 0.867 to 0.921. 
The mean estimate of 0. 894 with a standard devia- 
tion of 0. 009 for this range. The difference in k e 
for this shaft and the test specimens is negligible. 
All other modifying factors are assumed to be unity. 


*Normal distribution notation (mean, stand- 
ard deviation). 


Table 5 shows that Sg = (5.65x10 ; 3.2x10 
psi) for this material for r s = 1 . 03 and that 
there is little change for r g = «. 

Using the algebra of normal functions and 
Eq. (7), the modified endurance limit for this 
case can be calculated as follows: 

= (0.894; 0.009)(5.65X10 4 ; 3.2X10 3 ) 

C 

O 

S e = (5.08X10 4 ; 2.75X10 ) psi 

The bending moment, M, at the shoulder was 
estimated to be (126; 11.4) in. -lb. The bending 
stress, cr xa , a t the shoulder is calculated as 
follows 


xa 


_M__ (126; 11.4) 

L (0.098 d 3 ; 0. 000735 d) 
C 


_ / 1.28X10 3 . 116. 


d 3 


.6. i \ 
d 3 ) 


(8) 


The distance to the outer fibers C is given by 
(0.5 d; 0.0025 d). Standard deviation for this 
estimate is based on a manufacturer's guaranteed 
tolerance, Ad, of 0.015 d. The moment of in- 
ertia, I, is given by (0.049 d 4 ; 0.00295 d ) for 
I this configuration. The ratio I/C is given by 
(O'. 098 d 3 ; 0.000735 d 3 ). Substituting these val- 
ues yields the solution for bending stress given 
in Eq. (8). 

The shearing moment, T, transmitted by this 
shaft is given to be (1000; 80) lb-in. The torsional 
stress Txyjn at the shoulder is calculated as 
follows: 

t - T - __ apop,; 80 ) 

Xym 1 (0.196 d 3 ; 0.00147 d 3 ) 

C 



The von Mises-Hencky stress components 
are given as 
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s 


a 


' 1 . 28x1 (T 
^ d 3 " 


116. 4\ 
d 3 ) 


using Eqs. (2) and (8), and 


s 


m 


8. 83X10 
d 3 


3 



using Eqs. (3) and (9). 


r 


s 


s a 

— 0.144 
b m 


The above two stress distributions should be 
synthesized into the failure governing stress dis- 
tribution f(sf), along r s in Fig. 10. Along T s 



hence using the partial derivatives method (22, 
pp. 228-235), and the parameters for f(s a ) and 
f(s m ) given above, we get 



where 

9 Sf s a ^ 8s f s m 
3s a . sj Ss m Sf 

Consequently, 



Thus f(Sf) should be coupled with f(sf) to deter- 
mine the shaft diameter with the specified re- 
liability. Consequently f(Sf) should be determined 
next. A probabilistic Goodman diagram can be 
prepared from the research data. Fig. 10 is the j 
resultant diagram for this case. The top portion 
iof the Goodman surface is obtained by drawing 
three straight lines to represent a normal sur- ; ' 

face. The mean estimate loci is drawn as a solid i 

line - ordinate (5.08xl0 4 ; 0) and r g = 1.03 at | 

(5. 08xl0^; 5. 4xl0 4 ). * From here on data from j j 
the combined stress fatigue reliability research 
machines has not been reduced as yet. The mid- 
dle portion of the Goodman surface connects the 
unnotched ultimate strength (1 . 78x10 ; 2.5x10 ) 
to the r g = 1.03 points. The right hand boundary 
is defined by the surface boundary (1.71x10°; 1 

3xl0 3 ) and the angle (tan -1 2 = 63. 5°) as the ab- ‘ 
jscissa scale was chosen to be two times the or- j j 
dinate scale to improve accuracy. A second sur- ! 
face is constructed to define S a and S m for 
this case. This surface is defined by a point 
(0; 0) and three lines passing through the origin 
subtending an angle of (16. 1±1°) with the abscissa, j 

*This point is based on the experimental 
data given in Table 5. 

This is the currently accepted practice. 

It's validity shall be established at the conclusion 
of the present phase of this research. 
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Neglecting the effect of the very small variability 
in r g) the intersection of the strength surface 
with the r g line defines S a to be approximately 
(1.91 xlO 4 ; 1. 10x10^) and . Sjp to be approxi- 
mately (1 . 31xl0 6 ; 1.21X10 3 ). Combining these 
two distributions gives a failure governing strength 
distribution along defined as (1 . 33x10 J ; 
1.23x10 s ). 

The mean and standard deviation of f(Sf) and 
f (Sf) have thus been determined. If it is assumed 
that these distributions are normal, then the pa- 
rameters of these distributions are known and 
are those determined' previously. The amount by 
which these two normal distributions overlap has 
been specified by a part reliability requirement 
of 0.999. Eq. (6) may be used to approximate 
the reliability of this complex case. Then for 
R = 0.999, -f/o-£ should be 3. 09 (21, Table 3, 
p. 351) or since 


t ^Sf “ ^Sf 



( 10 ) 


substituting previously obtained values into Eq. (10) 
and solving for d yields 


! d 6 - 0.134284 d 3 + 0.00412 = 0 

Solving for d gives a d of 0.438 in. Use a 
shaft diameter of 15/32 in. It must be pointed out 
that the portion of the Goodman diagram in 
Fig. 10 used to determine f(Sf) is conservative i 
hence this shaft's weight may be reduced further. 

The results for the probabilistic design 
; method are tabulated in Table 7 . Three modifying 
! factors are not required as the experimental data 
includes these factors. An a priori reliabililty is 
designed into this part by making use of difference 
statistics. The probabilistic design methodology 
results in better than a 38 percent savings in ma- 
terials for this case. Furthermore the design 
; integrity is known to have a reliability of 0.999, 
such that on the average no more than one in a 
1000 such shafts will fail while functioning during 
their designed for mission. 


CONCLUDING REMARKS 

The Lewis Research Center in conjunction 
with The University of Arizona has undertaken 
the task of developing a mechanical design me- 
thodology to stati s' 'Tally match the available 
strength to the imposed stress. Analytical and 
experimental methods have been developed to 
reduce the number of modifying factors necessary 
to determine a material's endurance limit. Im- 
portant design variables are handled as statistical 
distributions. Provisions for designing an a 
priori reliability into a mechanical component 
are explained. In order to work this methodology 
out completely for one case, methods for deter- 
mining (1) statistical strength distributions for 
SAE 4340 steel; (2) statistical stress variations 
for a rotating shaft subject to reverse bending 
and torque; and (3) statistical methods for relating 
strength to stress are given in this paper. 

The sample mechanical design problems 
given in this paper show that weight savings are 
possible when the new design methodology is 
compared with the existing von Mises-Hencky- 
Goodman design procedure. 

There are some serious problems in trying 
to put the probabilistic design methodology into 
broad usage. The most serious of these problems 
is the gross lack of distributional data for ma- 
terial strengths. A second problem is the lack 
of loading distributions and combinations of stress 
distributions. Experimental test programs are 
currently under way to obtain more data. These 
efforts are directed toward alloy and carbon steels 
suitable for use in the design of machined parts 
subjected to fatigue. 

A third problem is working with the tails of 
statistical functions. For valid work to be ac- 
complished in this region, the statistical functions 
must be based on an unbiased sample of the pop- 
ulation. One of the important tasks defined in 
Phase II of this research is to determine the 
effect of sample size on the discrimination ability 
of selecting the best distribution from among 
various distributions that fit the experimental 
data. 
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Table 1 - Operational Specifications for the Strength Distribution Research Machines 


£ 

Test Machine 

Accomodate a specimen rotating at 1800 rpm; produce and hold a 

Loading Mechanism for 

steady torque and a reversed bending moment; holding chuck 
1.0 in. diameter maximum; simple design employing "off-the- 
shelf" components. 

Simple device to produce, hold, and transmit desired steady 

Steady Torque: 

torque of 5400 lb-in. to test specimen. 

Loading Mechanism for 

Simple device to produce a reversed bending moment of 3450 

Reversed Bending: 

in. -lb while specimen is rotating. 

Test Specimen: 

SAE 4340 steel, Condition C-4; MIL-S-5000B, certification of 

Instrumentation: 
Strain gages 

chemical and physical properties; uniform quality, same heat 
and processing, heat treat to Rockwell "C" 35/40 as per MIL- 
H-6875 with minimum tempering temperature of 1000° F; in- 
spection as per MIL-I-6868; D = 0. 735 in. , d = 0. 500 in. , r = 
0.150 in.; K b = 1.45; K s = 1.22, with 1 in. x 1/8 in. xl/8 in. 
keyways. 

To obtain dynamic arid static strain measurements in bending and 

Channels 

torsion. 

To handle at least 6 sets of strain gage outputs simultaneously. 

Slip rings 

To transfer strain gage data to amplifier while specimen is rotating. 

Amplifier 

To amplify static and dynamic output from strain gages. 

Recorder 

To produce a permanent record of amplified strain gage outputs. 


a Three such machines have been built to date and are being used to obtain distributional 
complex fatigue data for the steel specimens specified above. 




Table 2 - Cycles-To-Failure Fatigue Data Material: SAE 4340 Steel, Rockwell 35/40 


MIL-S-5000B, Condition C4 


r s< a ) 

CO 

0.825 

f**s 



s m = ° 



0.876 

0.800 

0.818 

0.805 

ffr s 






0.226 

0.082 

0.070 

0.056 

K V 

* 





3.1 

8.0 

6.2 

4. 3 

s a 

150,000 

120,000 

100,000 

80,000 

60,000 

120,000 

100,000 

80, 000 

75,000 

ps a (b) 

154.0 

114.0 

98.0 

81.5 

73.0 

111.0 

91.5 

75. 5 

65.0 

ns a (c) 

1.5 

0.0 

2.7 

0.9 

1.8 

1.3 

6.7 

3. 1 

3. 9 

Ky(d) 

1.0 

0.6 

2.8 

1.1 

2.5 

1.2 

7.3 

4.2 

6. 0 

n 

1 

1,550 

7,000 

16,250 

60,000 

93,300 

5,400 

12,950 

39, 550 

99, 400 

2 

1,950 

7,600 

16,500 

64,850 

103,400 

5,650 

13,050 

43,950 

100,900 

3 

2,450 

7,700 

16,900 

65,000 

124,150 

5,900 

13,350 

46,450 

101,200 

4 

2,750 

8,000 

18,150 

67,750 

125,950 

6,050 

13,850 

47, 550 

103,850 

5 

2,900 

8,100 

19,350 

71,150 

127,800 

6,100 

14,900 

52,700 

104, 350 

6 

2,950 

8,400 

20,600 

71,550 

145,800 

6,100 

16,700 

58, 350 

110, 850 

7 

2,950 

8,850 

21,100 

71,950 

155,200 

6,150 

19,100 

60,000 

117,900 

8 

3,050 

8,950 

21,200 

73,100 

172, 300 

2,200 

20,200 

60,000 

120,550 

9 

3,050 

9,100 

21,200 

74,700 

172,550 

6,900 

20,200 

61,800 

129, 250 

10 

3,200 

9,250 

22,550 

74,700 

172,650 

7,800 

20,450 

62,250 

131,950 

11 

3,200 

9,300 

22,900 

75, 650 

177,500 

7,800 

21,000 

65,850 

135,950 

12 

3,250 

9,350 

23,650 

78,850 

178,200 

8,600 

23,100 

67,500 

137, 350 

13 

( e ) 

9,750 

24,300 

82,250 

182,850 


24,350 

67,600 

140,350 

14 


9,800 

25,200 

83,800 

183,300 


24,650 

69,150 

142,400 

15 


10,000 

27,000 

86,350 

195,800 


25,000 

70,200 

143,200 

16 


10, 350 

27,250 

96,000 

196,900 


25,150 

72,000 

149,950 

17 


10,350 

27,450 

98,450 

203,000 


28,250 

75,800 

159, 300 

18 


10,550 

27,550 

107,400 

205,050 


32,100 

77,850 

167,750 


Additional stress ratios of 3. 0 and 0. 3 are being investigated. 

^Mean stresses rounded to the nearest 500 psi. 

^Variability cr sa is the standard deviation of the distribution about the mean. 
Coefficient of variation Ky is the ratio of the variability to the mean expressed as 
a percentage 

e Cycles-to-failure rounded to the nearest 50 cycles. 































Table 3 - Tensile Strength 
Data for Notched Specimens 3 


Table 4 - Tensile Strength Data for Unnotched Specimens 








Sut 

> 

o 

°b> 

Test 

Kut> 

Lb 

> 

kpsi 

kpsi 

No. 

klb 

klb 

(b, c 

) 

(b,c) 

1 

2 

ESmII 

47 

47 

0 

0 

253 

255 

5 

0 

305.0 

305.0 

3 

49. 

4 

46 

3 

254 

0 

299. 5 

4 

50. 

3 

47 

4 

359 

0 

299. 5 

5 

48. 

8 

46 

0 

251 

0 

306. 5 

6 

49. 

2 

36 

0 

253 

0 

302.5 


49. 

6 

46 

8 

255 

0 

304. 5 


49. 

8 

47 

1 

256 

0 

305. 5 

. 9 

50. 

5 

47 

7 

260 

0 

309.5 

10 

49. 

_9_ 

47 

A 

256 

_5_ 

302. 0 


a Specimen diameter at the base 
of the notched was controlled 
to 0. 4975 ±0 . 0025 in. which 
gives an area of 0.1944 
±0.0019 in. 2 . 

b Strength has been rounded to 
nearest 500 psi. 
c Data reduction to normal pa- 
rameters. 



*Yit> ®b> . 

kpsi kpsi 

(b,c) (b,c) 


31.5 32.5 24.3 0.4753 0.1774 177.5 183 

30.6 31.5 23.5 0.4764 0.1784 171 5 176 

30.6 31.6 23.3 0.4754 0.1775 172.5 178 

20.8 31.0 22.8 0.4755 0.1776 168.0 174 

29.9 31.1 22.9 0.4758 0.1778 168.0 175 

30.1 31.3 23.2 0.4722 0.1752 172.0 178 

30.4 32.5 25.0 0.4787 0.1800 169.0 181 

30.2 31.4 24.2 0.4757 0.1777 170.0 178 

30.6 31.6 23.6 0.4763 0.1782 172.0 177 

30.3 31.2 24.2 0.4755 0.1766 171.0 176 


a Specimen diameter was controlled to d ±0.0025 in. 
— o 

gives an area of A ±0.0019 in. . 

^Strengths rounded to the nearest 500 psi. 

c Data reduction to normal parameters: 

Yield Ultimate Breaking 

fi 171.0 178.0 255.0 kpsi 

a 3.0 2.5 4. 5 kpsi 


0 264.0 
5 254.0 
0 249.0 
5 251.5 
0 254.5 
5 256.5 
0 256.0 
0 253.0 
5 259.0 
5 250.5 


which 


Ultimate Breaking 
(i 255. 5 304.0 kpsi 

<j 2.5 3.0 kpsi 



































Table 5 - Endurance Strength Data 


Material: SAE 4340 Steel, Rockwell 35/40 
MIL-S-5000B. Condition C4 



OO 

Wp, 

R 

lb 

b = 1 

<b) 

c> = 7 

(a) 

25 

S 

26 

F 

25 

S 

26 

F 

25 

F 

24 

S 

25 

F 

24 

S 

25 

S 

26 

F 

25 

S 

26 

F 

25 

S 

26 

F 

25 

S 

26 

S 

27 

S 

28 

F 

27 

F 

26 

S 

27 

F 

26 

F 

25 

S 

26 

F 

25 

S 

26 

F 

25 

S 

26 

F 

25 

F 

24 

F 

23 

S 

24 

F 


Data 26 S Data Tensile testing 

Needed 27 F Needed no data required 


“Test results: 

S - success, if specimen survived bxi0 c ' cycles 
F - failure, if specimen broke before bxlo c cycles 


°Data reduction to normal parameters: 


r s 

CO 

3.0 

1.03 

0. 3 

0 


56.0 


56. 5 


Tensile 























Table 6 - Strength Distribution Best Fit Summary 

Material: SAE 4340 Steel, Rockwell 35/40 
MIL-S-5000B, Condition C4 


r s 

s a 

kpsi 

(a) 

Eqs. 

(b) 

r 

(cj 

Normal 

Lognormal 

Normal 

Lognormal 


144.0 

Y=490x+2780 " 

Y=6. 39X+8. 14 

0.9022 

0.8724 


114.0 

Y=1070x+9030 

Y=7 . 02X+9. 09 

0.9869 

0.9821 

OO 

98.0 

Y=3950X+22, 170 

Y=8. 37X+10.03 

0.9797 

0.9772 


81.5 

Y=1 2 , 7 50X+7 7,980 

Y=9. 50X+11.28 

0.9612 

0.9786 


73.0 

Y=34, 900X+161 , 980 

Y=10. 65X+12. 00 

0.9631 

0.9426 

3. o 

Values Needed 


111.0 

Y=1000x+6550 

Y=6. 95X+8 . 77 

0.9291 

0.9435 

fi ft 9^ 

91.7 

Y=57 80X+20 , 47 0 

Y=8. 79X+9. 92 

0.9779 

0.9752 


75.5 

Y=ll , 550X+61 , 030 

Y=9. 52X+H.03 

0.9818 

0.9665 


65.0 

Y=21 , 930X+127 , 580 

Y=10. lx+ll. 77 

0.9777 

0.9763 

0. 3 

Values Needed 

0 

0 

Y=255, 360X+2830 

Y=12. 48x+7. 97 

0.9891 

0.9894 


a Stress rounded to nearest 500 psi. 

■L 

°Equation of a straight line on special graph paper with ordinate (Y), 
abscissa (X), slope (m) and intercept (b) such that Y=mX+b. 
c The coefficient of correlation between two perfectly correlated linear 
variables should be 1.0000. Deviation from this value is a measure 
of goodness of fit. 































Table 7 - Alternator Inner Rotor 
Shaft Design Summary 


Parameter 

Conventional 

Probabilistic 

k a 

0.65 

N.R.^. 

k b 

0.85 

N.R. 

k c 

0.752 

D.M.O 3 ) 

k d 

0.867 

(0.894; 0.009)( a ) 

k e 

0. 680 min 

N.R. 

j Se 

9. 35xl0 4 psi 

(5.65X10 4 ; 3.2xl0 3 )psi 

: S e 

2. 28X10 4 psi 

(5.08X10 4 ; 2.75xl0 3 )psi 

: T 

1000 lb-in. 

(lxlO 3 : 80) lb-in. 

1 M 

126 in. -lb 

(1 . 26x10 ; 11.4) in. -lb 


1.28X10 3 j 

1.28X10 3 . 116.4 

b a 

d^ 

-3 ’ -3 P 

d 3 d 3 

s m 

i 

8. 8 3x10 3 pol - 
d 3 , 

8.83X10 3 . 709 

-3 ’ -3 P 

d 3 d 3 

; s a 

12. 3x10 3 psi 

(1.91X10 4 ; 1. 10x10 s ) psi 


§m 

85.0X10 3 psi 

(1.31X10 5 ; 1.21x10 s ) psi 


d 

19/32 in. 

15/32 in. ( 


R 

(c) 

0.999 ! 


a Normal distribution notation (mean, standard 

deviation). { 

k Legend: | 

(N.R. - not required as experimental data includes 
this factor. 

|D. M. - a different method is used to include this 
| factor. 

c The reliability obtained by the approach of Shigley 
|(14, p. 169) is not theoretically applicable to this 
(complex case, hence the reliability is really not j 
! known. 


* -j 


I 
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MATERIAL: SAE 4340 STEEL 

MIL SPECS: S-5000B, H-6875, 1-6868 

CONDITION C4, ROCKWELL C 35/40 



10 3 10 4 10 5 10 6 10 7 


CYCLES TO FAILURE, N f 

<B) STRESS RATIO: r s = 0.825 FOR CYCLES-TO- 
FAILURE AND r s - 1.03 FOR ENDURANCE. 

Figure 1. - Distributional S-N diagrams. 


r ALTERNATING BENDING STRENGTH DISTRIBUTION 

/ 

/ CONSTANTS: 



Figure 2. -Three-dimensional Goodman fatigue diagram. 
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y 



ff y * °ym * °ya 
T m " T )^m ± T xya 


Figure 3. - Rotating shaft surface element stresses. 



ds 

STRESS OR STRENGTH, kpsi 


Figure 4. - Probability areas of stress and strength. 
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Figure 5. - Stress and strength distributions. 



Figure 6. - A close-up view of the complex-fatigue research machines. 
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(A) DIMENSIONS, FINISH, FORCES AND STRESSES. 



(B) CONFIGURATION BEFORE AND AFTER FATIGUE TESTING. 
Figure 7. -Test specimen details. 



(A) INNER SHAFT DETAILS. 



(B) BENDING MOMENT DIAGRAM. 
Figure 8. - Alternator inner rotor shaft. 
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